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THE DETERMINATION OF DIPOLE REACTANCE BY 

PATTERN INTEGRATION 



RADIATION 



SUMMARY 

A new method for determining the input reactance of a cylindrical dipole 
is given based upon the principle of analytic continuation of Poynting integ- 
ration into the part of the radiation pattern outside the normal, or 'visible' 
region. The theoretical basis for this extension, recently described by 
Rhodes for the case of planar aerials, is re-stated here and then applied to 
cylindrical systems. The reactance of a thin dipole is evaluated for a range 
of dipole sizes and the results compared with those obtained by the well- 
known induced e.m.f. method. 



LIST OF PRINCIPAL SYMBOLS 



a dipole radius 

d separation of current filaments 
E electric intensity vector (volts/ 
metre) 
F(/x) Fourier transform of current distri- 
bution /(z) 
H magnetic intensity vector (amps/ 
metre) 
H^i\u), Hr(u) Hankel functions of the 2nd kind, 
zero and first orders (= Jo.i(w) - 
yYo.i(")) 
/ filamentary current 
/o dipole loop current 
/g uniform sheet current 
/(z) current distribution 
J^ current line density (amps/metre) 
k phase constant {= Irr/K) 
I dipole half-length 
P power (watts) 
P(a), P(/3) angular spectra 
P(V^), p(/^) pattern functions 



R 

u 

V 

W 

X, y, z 

X 

Z 

a, /3 

Hy) 
log 7 

Vo 

k 
- }k 



spherical co-ordinate variables 
radiation resistance (ohms) 



parameter defined by u = 

Ml - A^')'''^ 

parameter defined by v = 

ka{^^- 1)'/^ 

strip width 

cartesian co-ordinate variables 

input reactance (ohms) 

total input impedance ( = R + y'X) 

mutual and self-impedances 

angular spectrum variables 

Dirac delta function 

Euler's constant ( = 0*5772. . .) 

wave impedance of free space 

( 12077 ohms) 

wavelength 

pattern function variable (= sin a) 

cylindrical co-ordinate variables 

pattern function variable (= sin /3) 

scalar potential function 



1. INTRODUCTION 

There are two well-known methods for determin- 
ing the radiation resistance of an aerial. In the 
'Poynting vector method', integration of the far- 
field power pattern over all real angles gives the 
real (resistive) power radiated by the aerial. In the 
'induced-e.m.f. method', integration of the power 
flux over a surface tightly fitting the radiating ele- 
ments of the aerial gives the complex radiated 
power; the real part is identical to that given by 



the Poynting method whereas the imaginary part 
corresponds to the power stored in the so-called in- 
duction field of the aerial and this determines the 
input reactance. Since the Poynting vector method 
uses the far-field expressions for pattern integration 
it has always been assumed that it cannot be em- 
ployed to determine aerial reactance. In a recent 
paper by Rhodes, however, it was shown that by ex- 
tending the range of integration into a complex 
angular region, the method will give both the radia- 
tion resistance and the input reactance of a planar 



aerial. The theoretical basis for this extension of 
the classical Poynting vector method to the com- 
plex angular region, or 'invisible' part of the radia- 
tion pattern, had been set out previously in a most 
elegant paper by Woodward and Lawson . Their 
starting point was the plane-wave integral represent- 
ation^ of electromagnetic fields existing on one 
side of a plane aperture. It was shown that, for 
transverse-magnetic (TM) fields in two dimensions, 
the far-field radiation pattern is identical to the 
'angular spectrum' of plane-waves in the visible 
region (i.e. corresponding to the range of real an- 
gles). It was also shown that the radiation pattern 
and the fields at the aperture are related through 
the Fourier transform, and finally, that the radia- 
tion pattern lying outside the visible region deter- 
mines the reactive power stored in the vicinity of 
the aperture. 




Fig. 1 - Angular spectrum of plane waves 



to OX, with magnitude P(/3) d/3 referred to the 
origin, 0, is fully specified by the expression: 

P{B)d(ie "'^^^ '^°^^ * ^ sin^) 



In Reference 1, these results were extended to 
the general case in which the total fields of a lin- 
ear, time-invariant, planar aerial are represented by 
a combination of TE - and TM - fields. The analy- 
sis was taken one step further by integrating the 
Poynting vector over the visible and invisible parts 
of the radiation pattern to determijie the input im- 
pedance of a planar aerial. 

The main purpose of this report is to extend the 
new principle to cylindrical systems and, in partic- 
ular, to employ it for determining the reactance of a 
cylindrical dipole. As the method is not widely 
known, the report commences with a tutorial intro- 
duction to the basic principles involved. 



2. BASIC PRINCIPLES 

The theoretical basis for the extension of the 
classical Poynting vector method to the invisible 
part of the radiation pattern is most easily demon- 
strated for the case of planar aerials; the key to 
the method is the concept of an angular spectrum of 
waves and its identity with the far-field radiation 
pattern. 

2.1. Angular Spectrum and Radiation Pattern 

Following Woodward and Lawson , we con- 
sider the two-dimensional system shown in Fig. 1 
in which OYZ is the aperture plane with waves be- 
ing propagated only into the half-space given by 
X > 0. There is no field variation in the z-direction 
and the electric field E is chosen to be everywhere 
parallel to OZ; symmetry between E and H enables 
7]oW to be interchanged with E throughout, where 
7]o is the wave-impedance of free-space. 

An infinite plane-wave directed at an angle /3 



in which k = 27T/>vwith >v the free-space wavelength; 
the time-factor e'*"' has been suppressed*. The 
function P(/3) represents the angular spectrum of a 
set of infinite plane-waves and is a complex quan- 
tity. In general, the fields at a point (x, y) will be 
due to a combination of such waves so that we have: 



Ez{x,y) = I Pi/3) e 



/' 



-jkCx cos/3 + y sin/3) 



d/3 



(1) 



The magnetic intensities Hx(x,y) and Hy(x,y) are 
given by expressions similar to Equation (1) through 
the relations rjo dHx = sin/3 dEz and rjo dHy = 
- cos/3 dEz; these expressions are: 

Hx(x,y) =± I P(/3)e~'''^^ ^°^^ + ^ ^i"^> sin/3 d/3 
rjo.J 
c 

(2) 



and: 



fiy(x,y) =~— I P(/3)e "''"'''' ^°^^ ^y s'"^> cos/3 d/3 
Vo J 



(3) 



The remaining fields are given by Ex = Ey = Hz - 0. 

In the above equations, the contour of integra- 
tion C in the /3-plane must be carefully selected; 
it must be such that, for positive x and for positive 
and negative y, the exponential term does not tend 
to infinity at any point along it. We may therefore 
admit complex values of /3 {= /3r + j/3i) provided 
that cos/3 is given no positive imaginary part and 
that sin/3 remains purely real. It follows from these 
conditions that the signs of /3r and /3/ must be 

* Since the time-factor is e ' a positive travelling 
wave is given by e ' , in Reference 3, the negative 
exponent is used for the time-factor. 



identical (|/3r|< tt) and that /3/ may not be non- 
zero unless /3r = ± ~. The resultant contour of 
integration is shown in Fig. 2(a); path C2 is over 
the range of real angles (visible pattern region) and 
paths Ci, C3 are over the allowed range of complex 
angles (invisible pattern region). 

Equations (1), (2) and (3) give the field con- 
ditions at any point in the region to the right of the 
aperture (x positive) including the aperture plane 
itself (x = 0). By evaluating the field quantities at 
a great distance from the aperture we can show that 
the angular spectrum P(/3) is just the far-field com- 
plex radiation pattern of the source. Selecting 
Ez(x,y) and referring to Fig. 1 we can re-write 
Equation (1) as: 

Ez(r,0) = j P(/3)e-''" "^"^^"^ "^^ d/3 (4) 
c 
where r is the radial distance from the origin 0. For 
large r, this integral may be evaluated by tlie meth- 
od of stationary phase , since the integrand makes 
a significant contribution only near the region /3 = 
0; moreover, the contour paths Ci, C3 are excluded 
from the range of integration since the integrand is 
greatly attenuated in the invisible region with r 
large. Near /3 = 0, we may replace P(/3) with P(0) 
and expand cos(0 -/3) into 1 - (0 - /S)*/2; Equa- 
tion (4) then becomes: 



EziT,4>)^P{4>)e 



-ikr r 
-Tr/2 



i^ict> - /3) = 



d/3 



The limits of integration may be extended to + co 
without invalidating the result because of the high- 
ly oscillatory nature of the integrand. Evaluation of 
the resulting standard integral gives: 



Jh 



fCa (visible) 5- 
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^ C3 (invisible) 
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(b) 

Fig. 2 - Contours of integration 
(a) /3-plane ib) \p-p\ane 

and replacing P(./3) by -—■ P(./3) cos /3 in this equa- 

'70 . 

tion gives Hy(y). It is convenient to introduce a 

new variable defined by: 



= sin/3 



(7) 



Ez(.r,4>)^' 



'(r/A)" 



iq- - fcr) 



so that the path of integration in the /3-plane (Fig. 
(5) 2(a)) becomes the real axis in the (//-plane (Fig. 2 

(b)). Equation (6) can now be written as: 



which identifies P(0) as the far-field radiation patt- 
ern of the two-dimensional source. With a three- 
dimensional system, the above method yields the 
dependence for the far-field of a point 



familiar r ^ 



source 



1 



2.2. Fourier Relationships 

The fields existing at the aperture plane (Fig. 
1) are given by Equations (1), (2) and (3) with x = 0; 
i.e. the tangential electric iirjtensity is given by: 



E.iy) 



J 



P(/S)e 



-jky sin/3 



dyS 



(6) 



ez(y) 



CO 



pW 



(1 



■4>H^- 



i^y'f'dyij 



(8) 



in which p(>/') (=P(yS)) is a new pattern function 
(the angular spectrum expressed as a function of 
sin/3); the visible region now corresponds to 
-\<\p<\ and the invisible region lies outside this 
range, i.e. I'/'KI. We can now invert the relation- 
ship between £z(y) and p(>/') in Equation (8) by the 
Fourier transform theorem, viz. 



p{0) 



fcCl-V'V'^ 



2-n 



00 

f Ez(y)e>''y'^dy (9) 



Equations (8) and (9) are the connecting relations 
between the far-field pattern function p(0) and the 
tangential electric intensity £z(y) in the aperture. 
Starting from Equation (3), it can be shown that the 
corresponding pair of equations involving the tan- 
gential magnetic intensity in the aperture are: 



Hyiy) = 



CO 



p(0)e'''^>"^d0 (10) 



and 



CO 

p(0) = --^ f Hy(y)e'''y'^dy 



(11) 



The magnetic intensity Hyiy) is identical to J ziy) > 
where Jz(y) is the current density in a plane aeri- 
al*, which is also independent of z. Bearing in 
mind that, in practice, the width of the aperture in 
the y-direction would be finite we see that Equa- 
tion (10) can be used as the basis for synthesis of 
a continuous linear source from a given pattern 
function. 

We are now in a position to discuss the sig- 
nificance of the visible and invisible parts of the 
radiation pattern function. 

2.3. Poynting Integration with Complex Angles 

From Poynting' s theorem, the mean power flux 
density through any point in the aperture (Fig. 1) in 
the positive A;-direction is the product of -£z(y) and 
//y(y)**; hence the mean power P through each in- 
finite aperture strip of unit heighten the z-direction) 
is: 



P = 



00 

i 



Ez{y)Hy{y) dy 



brackets obtaining an integral expression for P in 
terms of the pattern function p(v/') (see Appendix, 
Section 7): 



277 

kri 



oj (1-^=)''' 



d4i 



(12) 



The factor 27rhas appeared in this equation because 
we have employed the asymmetrical transform (Sec- 
tion 2.2). Equation (12) gives the mean complex 
power flowing through each unit strip of the aper- 
ture; the real part Re(P) is the resistive (radiated) 
power and the imaginary part Im(P) is the reactive 
power, both quantities being expressed in terms of 
the far-field power pattern \p(\p)f. Because 
{l-xp"^)"- is real for -l<v/;<l and imaginary out- 
side this range, it follows that Re(P) is obtained by 
restricting the pattern integration to the visible 
region, -1<0<1, (the conventional Poynting vec- 
tor method), and also that Im(P) is obtained by 
allowing the integration to extend over the invis- 
ible region, |0|>1, where the denominator of 
the weighted power pattern, (1-0^)'-^, becomes 
-/(0'-l)'/M.e. 



Re(P) = 



277 r \pWr 

kvr, X a-V'^)'' 



Im(P) 



277 

kVo 



o -^i 



-1 



i^l 



d\p 






(13)* 



Once the pattern function |p(0) | of the strip is 
known. Equations (13) give the radiation resistance 
R and the input reactance X through the relations, 
Re(P) = I^R and Im(P) = I^X where Is is the curr- 
ent at the input terminals of the aerial occupying 
the aperture. Equations (13) can be combined into 
a single expression for the total input impedance 
Z(= R + jX) by writing 



Substituting the expressions for £z (y) and Hy(y) 
from Equations (8) and (10) we have: 



We can now apply the convolution theorem in order 
to invert the product of the transforms within the 



* We are only concerned here with the current which flows 
on one side of the "sheet". 

** The superior bar denotes the complex conjugate. 



277 r \pW\ 
kVols J (1-'/'^) 



dp 



(14) 



In using Equation (14) we recall tliat the denomin 
ator of the integrand becomes 
visible region, |0 | > 1. 



/(0^-l)'/^ in the in- 



2.4. Illustrative Examples 

Let us consider an aperture which is uniformly 
illuminated over a total width W (Fig. 3(a)), then we 
have Is/W = Jz(y) = tiy(y). In order to evaluate R 

* "Re" and "Im" denote the "real" and "imaginary part 
of respectively. 



and X it is first necessary to obtain p(V^) through 
the transform relation given in Equation (11) by sub- 
stituting Hy(y) = Is/W. The rectangular current dis- 
tribution transforms into a "sine" pattern function* 
giving, 



p(V^) = 



-Vols sinikW^/2) 



W^ 







z 










Ji.W 




1m 




y/Av//^/. 




_L_^ 


w 

2 







w 

2 


Y 



(a) 



The power pattern, |p(V^)r, illustrated in Fig. 3 
(b), can now be inserted into Equation (14) which 
yields the following integral expressions for R and 
X: 



R =■ 



4i1o r sin^(fcH'>/'/2) 



■nkW 



r sin (.kW^/2 
1 0^(1-0^/^ 



dp 



X = 



4i1o r si 
rrkW^ J 



sin (/cW 0/2) 



■nkW'^ J 0^(i//2-l)'/2 



■dv/; 



The above integrals, although quite formidable for 
arbitrary W, become straightforward for the extreme 
cases of an infinitely narrow (W-'O) and an infin- 
itely wide strip (H'-'Co) where the sine-squared pat- 
tern collapses into an omnidirectional and a delta- 
function pattern respectively. For example, in the 
first case with W-'O, we find that R = krio/2 and 
X = cx) so that a uniform filamentary aerial (radia- 
ting from both sides) has a radiation resistance 
dOrr/k ohms per unit length and infinite reactance; 
this result agrees with results obtained by conven- 
tional methods. 

As a second example, we shall show how the 
same technique can be used to determine the mutual 
impedance between two parallel radiating infinite 
filamentary currents / spaced d apart; the previous 
results maybe used directly if it is assumed initially 
that the filaments are plane and radiate from one 
side only. The magnetic intensity over the aper- 
ture may be written as, 

Hy(y) = l{h(y + d/2) +S(y-d/2)} 

where S(y) is the Dirac delta function, and if this 
is substituted into Equation (11) we obtain the pat- 
tern function p(i/'), i.e. 



p(V^) 



- _ gQgl 

V 2 , 



\pm' 




reactive power i radiatad power | reactive power 



(b) 



-e-iT -6fr -4ffi -2fr 
-1 



2fr i4fr 6fr Btr kW-^ 

1 



Fi^. 5 - Uniform current sheet carrying total current 

(a) current distribution J^(y) . 
(i>) power pattern function |p(>/')|^ 



The total complex power radiated by the aperture is 
2/ (Zg + Z{^) where Zs and Z;^ are the self-and 
mutual impedances of the filaments respectively. 
Equation (14) therefore gives an expression for 
(Zs + Zm) on inserting /| = 2/^ and \p(\p) |^ from 
the above relation, viz 



Zs +2m 



IkriQ r cos 



Hkd^/2) 
■ —dxP 



^y 



CO . 

krjQ f I + cos(/cdi//) 



02)' 



But we have already shown that for a current fila- 
ment. 



Zs 



vl 



d4: 



(l-02)'/2 



SO it follows that 



CD 

■1 



/cT)o ( cosikd^) 
ZM = I ™ dp 



^ \ (i-^y^ 



* The sine function is defined by the relation sine x 
sin rrx 

TTx 



The latter integral is in standard form and gives 



where H^ is the zero-order Hankel function of the 
second kind. We therefore have the result that the 
mutual impedance between two parallel current 
filaments (radiating from both sides) is: 



6^ ^(2) (^^ 



ohms per unit length which is in agreement with 
that obtained by conventional methods. 

It has been demonstrated that the total imped- 
ance Z of a two-dimensional plane aerial can be 
obtained by Poynting integration of the far-field pow- 
er pattern; the extension to the three-dimensional 
case is contained in Reference 1, where the resis- 
tance and reactance of a narrow planar dipole are 
evaluated by this method. 



3. APPLICATION TO CYLINDRICAL AERIALS 

The principles underlying the new method of de- 
termining the total impedance of a planar aerial from 
its radiation pattern were described in Section 2. 
We can now turn our attention to the problem of de- 
veloping a similar method for evaluating the reac- 
tance of cylindrical aerials. In this case, our start- 
ing point will involve the concept of an angular 
spectrum of cylindrical waves. The discussion 
which follows will be limited to systems possess- 
ing axial symmetry. 

3.1. Integral Relationships for Cylindrical Aerials 

We shall consider the three-dimensional system 
shown in Fig. 4 in which a cylindrical aperture of 
radius a and with its axis directed along OZ radi- 
ates outward travelling cylindrical waves. An ele- 
mentary TM-cylindrical wave with axial symmetry 
directed at an angle a to the horizontal with magni- 
tude P{a)da referred to the origin O is fully speci- 
fied by the following scalar potential function, 



-jkz sin o- 



d±{p,z) = P{a)daH^^\kp COS a.)e 



where the time-factor e''^^ has been suppressed. 
The elementary field components are then given 
through the expressions , 

ik 3 
dEz = k^ cos'^ ad^ , dH^= -— —-(d^) 

Vo dp ~ 

jk sin a o 

^g^ = + _____ ..__ (^0)^ ciEd, = dHp = dH^ = 0. 

p 'd P 

since there is no variation with 0. In general, the 
total field at a point (p,z) outside the aperture will 
be due to a combination of such waves so that we 



fz (A^ ) 



W<^(p,z) 




Fig. 4 - Angular spectrum of cylindrical waves with 
axial symmetry 

can write, 
Ez(P,z) = lp(a)k^ cos^aH^o^(kpcosa)e-i'^'^ ^i^^-da 



(15) 



and 

Hd> (A z) 



jk f 
— )P(a) 

"no -J 



)k. cosaHi\kpcosa)e-i'^^ ^^"°da 



(16) 



Equations (15) and (16) can alternatively be derived 
directly from Maxwell's equations in cylindrical co- 
ordinates. As with the plane-wave case, symmetry 
between E and H enables rjoH to be interchanged 
with E throughout so that the following analysis 
also includes systems involving transverse-electric 
(TE) cylindrical waves. 

The contour of integration C in the above in- 
tegral expressions must be carefully chosen; the 
path must be such that for p^a (positive) and 
. CD < z < CD the cosine-Hankel-exponential func- 
tions in the integrands of Equations (15) and (16) 
remain finite over the path in the a-plane. The con- 
dition involving the exponential function is as be- 
fore, namely complex values of a{= a/? + /a/) are 
admitted provided that sin a remains purely real, 
i.e. a/ must equal zero unless Of^ ~ ±tt/2; the 
conditions involving the Hankel functions are not 
so straightforward and will now be elaborated. The 



functions /i„ ,ikp cos a) must not increase indef- 

o » 1 ' 

initely as pi^a) is increased towards infinity for 
values of a which, from the previous condition are 
given by 



it may be shown from Equation (16) that 



VoH4>(r,6) = Eeir,0) 



a = ofi and a = ±— ± /a/; 

we need only consider, therefore, the argument 
kpcoso. in each of three states: finite real and 
positive, zero, and purely imaginary. The following 
properties of the Hankel function are relevant: 

(i) f^Q,J-^^ is complex unless v is imaginary when 
H^ iv) is imaginary but H^ (p) is real. 

(ii) For u equal to zero, \H^ ^(v)\ becomes infin- 
ite. 



so it follows that PW) is also the far-field radiation 
pattern for /i</,. 

There is a Fourier relationship between the radi- 
ation pattern and the fields existing at the aperture 
ip - a); they are not so straightforward as those 
given for the plane-wave case in Section 2.2 because 
of the presence of the Hankel functions in the cyl- 
indrical field expressions. It is convenient to in- 
troduce a new variable fd. defined by 



P- 



sina 



(18) 




(iii) The asymptotic formulae for large v are valid 
for V complex, i.e. 



"y'(i^ ~mT/ 2 -77/4) 



It follows that cosa must not contain any positive 
imaginary part, i.e. the signs of a/? and a/ must be 
identical (la^l^ tt); the poles in the Hankel func- 
tions which occur with a = ±tt/1 appear in the 
scalar wave function ^ but in the expressions for 
Ez and /i</> (Equations (15) and (16) they are anni- 
hilated by the cosa weighting functions. For the 
field quantities therefore, the resultant contour of 
integration in the Orplane is identical to that given 
previously for /S in the plane-wave case (Fig. 2 
(a)); as before, path C2 is over real (visible) angles 
and paths C\, C3 are over complex (invisible) an- 
gles. 

Evaluating the field quantities given by Equa- 
tions (15) and (16) at a great distance from the ap- 
erture (r-»oo) we can show that the angular spectrum 
P(a) is just the far-field complex radiation pattern 
of the source. Selecting E2,(p,z) and referring to 
Fig. 4 we can re-write Equation (15) in a form suit- 
able for stationary-phase integration (Section 2.1) 
by inserting p = rcosO, z = rsin6 and writing 
H^ ikpcosa) in asymptotic form for large r. On 
performing the integration we obtain the result. 



Ezir,e) =^2jk CO s e PiG) £;''"' 



so that the path of integration in the a-plane be- 
comes the real axis in the /Lt -plane (c. f. Equation 
(7), Fig. 2(b)). Equation (16), with p = a, can now 
be written as. 



H^(a,z) 



00 

^0 J 



p(fi}H{''hlcai]7J?)e''''^^dfM 



(19) 



in which p(a) (= P(a)) is a new pattern function; 
the visible region corresponds to -1< /Lt < 1 and the 
invisible region to |/l(,|>1. The relationship be- 
tween II^(.a,z) and pC/J,) in Equation (19) can be 
inverted by using the Fourier transform theorem 
which gives. 



00 



(20) 



Equations (19) and (20) are the connecting relations 
between the far-field pattern function p(/2) and the 
tangential magnetic intensity in the aperture 
//<^(a,z). There is a similar pair of relations in- 
volving the tangential electric intensity in the aper- 
ture Ezia,z) which can be derived from Equation 
(15); these relations are. 



(17) 



i 



Ez(a,z) = kH p(;Lt)(l-;Lt^)'/^M2)(^^Jl_^2)g-;fczM 



d/J. 
(21) 



which identifies PW) as the vertical radiation patt- 
ern of the source, i.e. 



Ee(.r,6) is proportional to Pid). Similarly, 



-1/ °° 

p(fL)= ^^''''^ ^ f Ezia,z)e'''^^dz 



(22) 
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If the cylindrical aperture of radius a, shown in 
Fig. 4, is replaced by an equivalent cylindrical 
sheet of current with distribution density J^^z), we 
cap use Equations (19) and (20) by writing H^(.a, z) = 
Jz(z). It follows that, if /(z) is the total current 
distribution on an axially symmetric cylinder, the 
equivalent tangential magnetic intensity at the sur- 
face is given by 

substituting this expression into Equation (20) gives 
the radiation pattern function of the cylinder. 

3.2. Poynting Integration 

The mean power flux density out of any point 
i4\z) in the cylindrical aperture (Fig. 4) is 



-Ez(a,z)H^(a,z) 

so that the total mean power P flowing outwards 
from the entire length of the aperture is: 



P= -2 



CO 

"1 



Ez^a,z)H^{a,z} dz 



Substituting the expressions for E^ and H^ from 
Equations (19) and (21) and applying the convolution 
theorem (c.f. Section 2.3) we obtain an integral ex- 
pression for P in terms of the pattern function P(yLi), 
viz 



Vo J 



pifi) 



H<^^kkah-fi'')H[^\ka h- u^) 



(l-u=)-'/^ 



d/J. 
(23) 



As a result of the contour of integration, the terms 
(1 -/i^)^ in Equation (23) must be taken as 

-yV=-l)'/^ 
in the invisible region. 

Finally, we can replace pi/x) in Equation (23) 
by the expression given in Equation (20) and obtain. 



CO 

vok r I 

iFCa) 
47T^a J 






d/J. 



(24) 



where F(fx) is the Fourier transform of the total 
cunent distribution /(z) along the equivalent cylin- 
drical aerial, i.e. 



F(/j) 



CO 

J 



,/ \ jkzu j„ 
/(z)e' dz 



(25) 



Equations (24) and (25) enable us to calculate 
the total power P radiated by a cylindrical aerial 
from its far-field radiation pattern; the weighted 
power pattern is integrated over the range of visible 
(- 1 < yLi < 1) and invi sible (\fx | > 1) angles. The res- 
istive (radiated) power is given by Re(P) and the 
reactive power by Im (P). In the visible region, 
both Hankel function expressions in Equation (24) 
are complex so that this region contains both radi- 
ated and reactive powers; in the invisible region, 
the zero-order Hankel function is imaginary but the 
first-order Hankel function is real so that this region 
contains only reactive power. This differs from the 
plane-wave case in which reactive power is deter- 
mined by performing the integration over the invis- 
ible region only. 

Equations (24) and (25) provide the basis for 
determining the total impedance of a cylindrical aer- 
ial of arbitrary length 21, radius a, and current dis- 
tribution /(z), provided that the integration can be 
performed. 

3.3. Thin-Cylinder Approximation 

If the radius a of the cylindrical aperture, and 
hence also of the equivalent cylindrical aerial, is 
assumed to be small, it is found that the integral 
in Equation (24) becomes more manageable. To 
show this we write the W^ and H^^' functions in 
Equation (24) as first order approximations for small 
(complex) arguments, i.e. 



o 
u -0 



Hy^'(u)^i~j^[og{m 



and 



/(2) 



'l\" (") - ik 



u -0 

where u is identical to ka(l-fx^)'^ in the visible 
region and -jkai/u.' -l)'-^ in the invisible region and 
log y(= 0*5772...) is Euler's constant. This step 
requires some justification; the visible and invis- 
ible parts of the integral must be considered sep- 
arately. 

In the visible region, u is real with a maximum 
value ka at /J. = 0. If ka = 0*3, the above approx- 
imations give maximum errors of approximately 
(2+y4)% in H^^^hu) and -10% in H[^hu) \ with 

\Y^iu)\>10J^(u). 



In this region therefore, it is reasonable to assume 
that the approximations are valid for ka < 0-3 . It 
follows that for cylinder radii a 4 0*05\, approximate 
integration of Equation (24) will yield sufficiently 
accurate results for the radiated power and also 
for that part of the reactive power which resides in 
the visible region. 

In the invisible region u becomes negative 
imaginary. Writing u = -Jv, the value of u ranges 
from zero up through large values, where vc^kafj. 
with fj. >> 1, to infinity. The above Hankel function 
approximations for small arguments now become, 



H(^>(-;v)^-;llog(^^ 



and 



H 



1 ' 7TV 

v-0 



In this case, v = ka {or fi ~ -J2 and if ka = 0-3 the 
above approximations give errors of approximately 
-/4% and -4% in H^*^(u) and H^f* {u) respectively 
at the invisible angle corresponding to /x = ^2; as 
yu, is increased above this figure, the approximations 
will become decreasingly accurate. However, the 
Hankel function expression in the integrand of 
Equation (24) is weighted by the factor v|F(yu,)p; 
it will be recalled that F(/x) is the Fourier trans- 
form of the current distribution along the cylinder. 
The rate at which this weighting factor falls off 
with increasing /x will therefore depend on the 
shape and extent of the current distribution /(z); 
narrow distributions give radiation patterns which 
extend far into the invisible region whereas wide 
distributions give relatively sharp patterns. We 
find that with cylinders of finite (but not short) 
lengths carrying sinusoidal current distributions, 
the transform F(yu,) is sufficiently sharp to validate 
the inclusion of the Hankel function approximations 
in the invisible part of the power integral for small 
cylinder radii. In such cases the tentative ass- 
umption* will be made that the approximate integrand 
will give sufficiently accurate results for total re- 
active power providing a-4:Q'Q5K. 

If the above approximations are made, the real 
and imaginary parts of the power integral (Equation 
(24)) can be written as, 

.« r 



and 



Equations (26) and (27) determine the resistance 
and reactance respectively of a thin cylinder of 
radius a from the Fourier transform F(/x) of its curr- 
ent distribution /(z) (Equation (25)). 

Alternatively, Equations (26) and (27) can be 
written in terms of the radiation pattern P(a). We 
note that the pattern function F(yu,)(l -yu,*) ''^ is prop- 
ortional to P(a) since yu, = sin a. It then follows 
that the result for radiated power is. 



Re(P) : 



n/l 



P(a) rcos a-da. 



~v/2 

where /f is a constant; this expression shows that 
Equation (26) is the conventional power integral. 

It will be seen that, with the thin-cylinder 
approximation, the resistive power is independent 
of cylinder thickness whereas the reactive power 
depends on thickness and length through a logar- 
ithmic factor in the integrand (Equation 27). 

4. REACTANCE OF A THIN DIPOLE 

The extension of Poynting integration into the 
invisible region has been developed for the case 
of cylindrical apertures; we now employ the method 
to determine the reactance of a thin dipole with a 
sinusoidal current distribution along its length. 
The result will be compared with that obtained by 
two other methods. 

4.1. Determination of Input Impedance 

We consider a thin dipole of radius a and 
length 11 with a sinusoidal current distribution 
i(z) as shown in Fig. 5. 



z-t 



Fig. 5 - 

Thin dipole with sinusoidal 
current distribution 



z--( 



/I(z)e^"' 



To 



2a 



Im(P) 



.^[iwiv-i){4...i.'-ik.o.Qj.. 



(27) 



* A rigorous discussion of the asymptotic behaviour of 
the integrand in Equation (24) is outside the scope of 
this report. 
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If we insert /(z) = /q sin fc(i - \z\) into Equation 
(25) we obtain its Fourier transform F(.fi), namely: 



Fi/J-) 



2/n 



fe(l-u=) 



(cos klfi - cos kl) 



(28) 



where /q is the loop current. The input impedance 
of the dipole Z(= R + jX) can now be evaluated by 
inserting Equation (28) into Equations (26) and (27) 
and using the relations, 






(29) 



Writing Tfc = 12077, the expressions for R and X 

become: 

1, 

cos kl)^ 
——d/j. (30) 

Vi - /i ) 





R on 120 cosec^fci 



1, 
j (cos kin • 

J (1 - 



and 



X ^' 



240 cosec^ k 
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Fig. 6 - Input reactance of a thin cylindrical dipole 

of radius a and length 21 determined by pattern 

integration 



cos kl/i - cos kl) 11 12 I 

— ■^— — — { - log yu. -1 + log 
(/x - 1) I 2 




(31) 



The integral in Equation (30) can be evaluated by 
standard methods and the result is, 



Rc± 



30 



sin^kl 



{2Ki(2/:0 + [Si(4/cO - 2Si(2fcO] sin 2kl 

- [Ki(4A:0 - 2Ki(2/c/)]cos 2kl} 

(32)* 



which is seen to be the usual expression for the 
radiation resistance of a filamentary dipole"^ of 
length 21; there is no dependence on the dipole 
radius a. 

The integral in Equation (31) is more difficult 
because the logarithmic part (first term) cannot be 
reduced to a standard form; the second term inte- 
grates into 77/4 log(7fea/2)sin(2fe/) so that the 
result is. 



60 i 2 f ( 
ain^kl J ttJ 



cos kl n - COS kl)^ 

(/x2 - 1) 



which is almost the same expression as that ob- 
tained for the narrow planar dipole of total width 
a; the two expressions differ only in that the latter 
result contains an additional term -180 cotkl which 
depends only on dipole length. For a half-wave 
dipole (kl = 77/2) the second temi in Equation (33) 
disappears and the reactance, as determined by 
numerical evaluation of the remaining integral, 
becomes 42'S ohms; this is in agreement with the 
value of 30Si(277) obtained for a cylindrical half- 
wave dipole by the induced e.m.f. method. 

For lengths other than odd multiples of a half- 
wavelength the input reactance depends on the 
radius as well as the length of the dipole. In Fig. 6, 
the reactance X determined from Equation (33) is 
plotted as a function of dipole half-length I for 
varying dipole radii a/k. The results show that the 
input reactance depends very largely on the dipole 
thickness and becomes very high for thin dipoles. 
The general variation is very similar to that ob- 
tained by the conventional induced e.m.f. method. 



[0g\j-L^- l|dyLi + log 



|Lf I sin(2fc0 



(33) 



Si(x) 



J sin V I < 

—dv and Ki(jc) = j 



{1 - cos v) 



-dv 
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4.2. Comparison with Other Methods 

It is of interest to compare the results obtained 
for the input reactance of a thin cylindrical dipole 
by pattern integration with those given by some 
other methods. Particularly important is the induced 
e.m.f. method which employs the result for the 
mutual impedance between two identical dipoles to 
determine their self-reactance. 
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Fig. 7 - Input reactance of thin cylindrical dipole 
comparison of results obtained from 
(a) pattern integration (b) induced e.m.f. method 



In Fig. 7 we have re-plotted two of the results 
from Fig. 6 together with those obtained by the 
conventional induced e.m.f. method* for the same 
dipole thicknesses. It will be seen that there is 
little difference between these results in the region 
of dipole-resonance {l/k ^ 0'25) but that for short 
dipoles the results diverge, with the present results 
(curves (a)) giving larger capacitive reactances 
than the induced e.m.f. method (curves (b)). (It will 
be recalled from Section 3.3 that the approximate 
theory (Equations (26) and (27)) is expected to give 
errors for short distributions.) For extremely thin 
dipoles, however, both solutions approximate to the 
expression, X ^ 120 cot{k[)[og{a/k), and the frac- 
tional difference between the results obtained by 
the two methods would be insignificant. 

The reactance of a thin cylindrical dipole may 
also be estimated from that of the planar dipole by 
using an electrostatic transformation which can be 
applied to uniform systems of conductors carrying 
TEM waves. Assuming that a cylindrical dipole of 
radius a is equivalent to a strip-dipole of total 



* The inducedfield is evaluated over a cylindrical sur- 
face of radius ^^2a (see Reference 8). 



width 4a, we find from Reference 1 that the new 
expression for input reactance differs from that 
given in Equation (33) by an additional capacitive 
reactance of approximately 14cot kl. This means 
that the difference between the two results is insig- 
nificant over the range of interest (i.e. near dipole- 
resonance) with the results diverging for very short 
dipoles; the latter effect is not surprising because 
the above transformation takes no account of end 
effects and the approximate power integration be- 
comes inaccurate for short distributions. 



5. CONCLUSIONS 

It has been shown that it is possible to deter- 
mine the reactance of a cylindrical dipole from its 
far-field radiation pattern by suitably integrating 
the power pattern over the visible and invisible 
angular regions; this has extended to cylindrical 
systems the new principle of analytic continuation 
of the Poynting vector method recently described 
for planar radiators. 

For thin dipoles, the Poynting integrand re- 
duces to a manageable expression and the reactance 
in this case has been evaluated for several dipole 
sizes and compared with that obtained using the 
well-known induced e.m.f. method. The conclusion 
is that the difference between the results of the two 
methods is insignificant over a large range of 
dipole sizes near the first resonance; the agreement 
is exact for the thin half-wave dipole and is inde- 
pendent of thickness. 

It is emphasized that the new method is not 
restricted to thin dipoles; the only difficulty in 
applying it to thick dipoles with arbitrary current 
distributions is actually evaluating the complex 
integral expressions for reactance which obtain in 
such cases. 

Finally, it should be pointed out that the new 
method, which involves complex power integration 
over the aerial surface, actually uses an induced 
e.m.f. technique; the essential difference between 
the new and the conventional methods is that in the 
former, mathematical manipulation allows the power 
integral to be stated in terms of the far-field radia- 
tion pattern of the aerial rather than its near-field 
expressions. 
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7. APPENDIX 
Collapse of Triple Power Integral into Single Form 



The triple integral to be considered is 



It therefore follows that 

CO 



/ = 



[| f/i(0)e"'''>''^#¥ f ^e'''>''^didy (34) J f,(4j)f^^p^)d^ = r^^T/, (0)7^0)} (36) 



where, from Section 2.3, we have 

/i(0) = (1 -^J^r'^'pi^) and Ui^P) = p(0) 
Equation (34) in its simplest form is written as: 



/ = 



00 



(.4i)TUi4i)dy 



(35) 



when T is an operator denoting "the Fourier trans- 
form of. The convolution theorem is now applied; 
in its usual form it can be written, 

00 

T I /i(0)/2(^ - 0) #= r/i(0) r/2(0) 



so that we have* 



00- 



^ /i (0)4(0 - x) d0 = r/i(0) r/2(0) 



* If r/(0) = Fis) and r/(0) = Gis), then Fis) = a~s). 



The r.h. s. of Equation (36) is the integral: 

00 

^ J r/i(0)r/2(0)e'''>'^dOcy) 

- 00 

If we now set x ~ 0, Equations (35) and (36) give 

00 

■^= J A(0)/2(0)d0 



so that we have. 



/ = 



00 



(4j)U(^P)d^ 



(37) 



Equation (37) is the single integral equivalent of 
Equation (34) and is tlie required result. This form 
of the power integral has been used in deriving 
Equation (12), Section 2.3, and Equation (23) 
Section 3.2 in the main text. 
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